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Abstract 

We conduct a rigorous investigation into the thermodynamic insta- 
bility of ideal Bose gas confined in a cubic box, without assuming 
thermodynamic limit nor continuous approximation. Based on the ex- 
act expression of canonical partition function, we perform numerical 
computations up to the number of particles one million. We report 
that if the number of particles is equal to or greater than a certain 
critical value, which turns out to be 7616, the ideal Bose gas subject 
to Dirichlet boundary condition reveals a thermodynamic instability. 
Accordingly we demonstrate - for the first time - that, a system con- 
sisting of finite number of particles can exhibit a discontinuous phase 
transition featuring a genuine mathematical singularity, provided we 
keep not volume but pressure constant. The specific number, 7616 can 
be regarded as a characteristic number of 'cube' that is the geometric 
shape of the box. 
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1 Introduction 



By definition, first-order phase transitions in thermodynamics feature a gen- 
uine mathematical singularity. Whether finite systems in Nature can liter- 
ally exhibit such an infinity is a long standing controversial question [Il|2]. 

When a thermodynamic system is composed of a definite number of 
particles, say A^, and is in contact with a heat reservoir, the key quantity is 
the canonical partition function: 

ZN{li,V) = Y.^e-P^^, (1) 

where the sum is over all the quantum states of the A^-body system. When 
the energy eigenvalues depend on the volume of the system, Zn is a func- 
tion of the volume V, and the temperature through j3 = l/{k^T), where 
/cb denotes Boltzmann constant. Provided the precise canonical partition 
function, we may compute various physical quantities, which include the 
pressure, the entropy, the internal energy and the specific heat per particle 
at constant volume as follows: 



P = {l/(3)dvlnZN, S = k^il-l3dp)lnZN, 

(2) 

E = -d^lnZN, Cv = ik^/N)(3^djlnZN . 

In particular, d'^ln Zn = {{E^ — E)"^) being a standard deviation squared, 
Cy is finite and never negative. Further, the temperature derivative of the 
probability for the system to occupy a certain quantum state ijj reads 



d_ 

dT 



V 



(e-P^^/ZN) = A:b/32 {E^ - E) e'^^^ jZ^ . (3) 



As the temperature increases from absolute zero to infinity, the correspond- 
ing probability increases if the energy eigenvalue is greater than the average 
i.e. E^ > E, and it starts to decrease in the opposite case, E^ < E. 

Since the canonical partition function appears as a positive definite an- 
alytical function of its arguments, clearly all the physical quantities listed 
in ([2]) and ([3]) should not feature any singularities. They may do so only in 
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the thermodynamic hmit: the limit of — t- oo and V ^ oo with N/V held 
fixed [3]. 

However, strictly speaking, infinite limits are hardly realistic and exist 
only in theory [Il[2l[31[5] . The above analysis seems to suggest that Nature 
does not admit a discontinuous phase transition featuring a genuine math- 
ematical singularity, which is somewhat different from experiments or our 
daily experiences under the standard pressure 1 atm. 

In this paper we pay attention to, among others, the fact that the above 
finiteness and continuity are for the cases of keeping the volume fixed. Once 
we switch to an alternative constraint of keeping the pressure constant, we 
demonstrate, for the first time, that canonical ensembles with finite number 
of physical degrees may undergo a discontinuous phase transition. 

The organization of the present paper is as follows: 
Section [2] is devoted to a systematic analysis on the thermodynamic insta- 
bility of a generic finite system. We explain how canonical ensembles with 
finite number of physical degrees may exhibit a discontinuous phase transi- 
tion when we keep not volume but pressure constant. Further we present 
a theorem which states that, a thermodynamic system must be unstable at 
low temperature near absolute zero if the ground state energy is volume 
independent. Examples include systems with vanishing ground state en- 
ergy, such as supersymmetric models, ideal Bose or Boltzmann gases under 
periodic or Neumann boundary conditions. 

In section [3l as a concrete model we focus on ideal Bose gas which 
is confined in a box and subject to Dirichlet boundary condition. Since 
there is no zero mode, the ground state energy depends on the volume 
and our theorem is not applicable in this case. Nevertheless, by numerical 
analysis we show that the ideal Bose gas reveals a thermodynamic instability 
and consequently undergoes a first-order phase transition, if the number of 
particles is equal to or greater than 7616. 

The final section [J] conveys our discussion. In particular, we comment on 
the similarity between the permutation symmetry of the identical particle 
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indices and the gauge symmetry in high energy physics. 

Appendix contains our numerical verification that ideal Bose or Boltz- 
mann gases under periodic or Neumann boundary conditions exhibit a ther- 
modynamic instability at low temperature near absolute zero. 

Although ideal Bose gas has been studied for decades [TBHlT] and dis- 
cussed in many textbooks [6IIT8H20] , the implication of the constant pressure 
constraint to the canonical ensemble of finite has been rarely explored^ 
To the best of our knowledge the finite A-effect on the canonical ensemble 
has been addressed only in the case of keeping the volume fixed, rather re- 
cently by Kleinert |23] and by Glaum, Kleinert and Pelster j24] . The earlier 
focus was typically on either grand canonical or micro-canonical ensembles, 
and the computations often assumed a continuous approximation to convert 
discrete sums to integrals |25p26j , unless an external harmonic potential sets 
the sum to be taken over a geometric series [271436] . 



^The so-called constant pressure ensemble [21II22| fixes the pressure and allows the 
volume to fluctuate, as its partition function is given by 

Yn{P,P)^ J dVe-^''''ZNW,V). 

However, just like Cv = {ke / N) 13'^ In Zn, the specific heat therein is positive definite 
and finite: 

< Cp = {ks/N)0^dl IuYn <oo. 

As we are interested in the precise change of the volume without allowing any fluctuation 
of it, in the present paper we focus on the standard canonical ensemble. 
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2 First-order phase transitions in finite systems 



Here we explain how canonical ensembles with finite number of physical 
degrees may exhibit a discontinuous phase transition featuring a genuine 
mathematical singularity when we keep not volume but pressure fixed. 

Prior to rigorous analysis, a thought experiment may provide an intuitive 
clue for this claim: If we fill a rigid box with water to the full capacity and 
heat it, the temperature will increase but hardly it evaporates. However, 
once we open the lid, simply it boils. 

Explicitly, as pressure being a function of T and V, from ([2]) we have 

dP = dTdv S + dVdv P. (4) 

Hence under constant pressure the temperature derivative acting on any 
function of (3 and V can be computed as 



d_ 

dT 



= -k^p'dp - (dvS/dvP) dv . (5) 
Note that unless explicitly specified as ^ and ^ , throughout the 



V 



paper dp always denotes the beta derivative at fixed V and dy is the volume 
derivative at fixed /3, as already taken in ([2]). 

In particular, the specific heat per particle under constant pressure reads 



r — Xd_ 

~ N dT 



^{E + PV) = Cv-{dvSy/{k^f3NdvP). (6) 



We pay attention to the denominator here which is essentially the volume 
derivative of the pressure: 

dvP = I3{ {dvE^ - {dvE^yf ) - {dlE^) . (7) 

This quantity possesses an indefinite sign. If dyP < 0, the system is stable: 
it resists against the change of external pressure by adjusting its volume. 
On the other hand, the opposite case, dyP > 0, characterizes the first- 
order phase transition, as the volume at the phase transition is not single 
valued [6]. Clearly from ([5]), when dyP = 0, a singularity develops and 
every physical quantity should change discontinuously. Moreover, when dyP 
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crosses the vanishing hne, Cp can be negative. To the best of our knowledge, 
all the known systems revealing negative specific heat do not keep the volume 
constant, such as in astrophysics [71[8], melting transitions [9|fT0]. and in real 
experiments pTt[T2]. 

We emphasize that from ([5|) the only way for a finite canonical ensemble 
under constant pressure to reveal singularities is through dyP = i.e. the 
sign of the thermodynamic instability. 

We close this section by presenting a sufficient, yet unnecessary, condi- 
tion for the thermodynamic instability. 



• Theorem 

At low temperature near absolute zero, a thermodynamic system is un- 
stable i.e. dyP > 0, if the ground state energy is volume independent. 

The proof is simple once we spell the canonical partition function as 
Zj\[ = J2n^ne~^^" ^ where E^s are the possible energy eigenvalues and Qn 
is the corresponding degeneracy. By direct manipulation we may express 
dyP in the following form: 

(8) 

where Eq is the volume independent ground state energy, dyEo = 0, with 
the degeneracy ^Iq; Ei is the first excited state energy satisfying dyEi / 
with the degeneracy Qi. The ellipsis denotes exponentially suppressed terms 
for large (3 or low temperature. Clearly at low temperature near absolute 
zero dyP becomes positive. This completes our proof. 



Examples include systems with vanishing ground state energy, such as 
supersymmetric models, ideal Bose or Boltzmann gases under periodic or 
Neumann boundary conditions. Appendix contains a numerical verification 
of the latter. 
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3 Ideal Bose gas confined in a cubic box 



Based on the general analysis of the previous section, henceforth as a con- 
crete model we focus on ideal Bose gas confined in a cubic box and subject 
to Dirichlet boundary condition. Since there is no zero mode, the ground 
state energy depends on the volume and our theorem above is not applicable 
in this case. Nevertheless, by numerical analysis we show that dyP therein 
assumes positive values for some interval of temperature if > 7616. 

3.1 Algebraic analysis 

For non- interacting identical bosonic particle systems, more easily computed 
than the canonical partition function is the grand canonical partition func- 
tion: 

= Un{l-r]e->'''^y' , (9) 

where rj denotes a fugacity and n corresponds to a good quantum number 
valued "vector" which uniquely specifies every quantum state of the single 
particle system. Taking logarithm and exponentiating back, we acquire an 
alternative useful expression: 

2 = exp(Er=i Afcr?Vfc) , A, := En e-'^^'^^ . (10) 

From the power series expansion of this, Z = Eat Zj^rj^ , one can easily 
read off the canonical partition function, as previously obtained by Matsub- 
ara [37] and Feynman [38J : 

Zn = Y.m. nil {Kr^/{ma\ a^^) , (11) 

where the sum is over all the partitions of A^, given by non-negative integers 
nia, a = 1,2, - ■ ■ ,N satisfying N = E^i a™a- In particular, with Ai = Zi the 
partition as mi = N leads to a conventional approximation |39], or the canonical 
partition function of ideal Boltzmann gas: 

ZNO^{Zif/m. (12) 

This approximation would be only valid if all the particles occupied distinct states, 
as in the case of high temperature limit. Other partitions then give corrections to 



2 = Un ET=oive 
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such underestimation: Compared to ideal Boltzmann gas, ideal Bose gas has higher 
probability for the particles to occupy the same quantum state. 

Yet, according to the Hardy-Ramanujan's estimation, the number of possible 
partitions grows exponentially like e'^ v^^/(4a/3A^), and this would make any 
numerical computation practically hard for large N. Alternatively, we consider a 
recurrence relation on Zn, which was first derived by Landsberg 1401 and can be 
easily reproduced here after differentiating Z in (fTOl) by 77 : With Zq = 1 we get for 
N>1, 

Zn = (Ek^iXkZN-k) /N. (13) 

Further, if we formally define an 00 x 00 triangularized matrix A whose entries are 
given by 

{AQ-b/(a-l) for a>b 
(14) 
otherwise , 

the above recurrence relation gets simplified: 

Zn = EZo A[iV+l, n+1] Z„ , (15) 
such that it has a solution given by a particular entry of a certain matrix: 

zn = (Er=o ^"m+i, 1] = (/ - A)-i[iv+i, 1] . (16) 

Since A'''[A^+1, 1] vanishes for k > A^+1, the triangularized matrix A can be effec- 
tively - and happily - truncated to its upper left (7V+l)x (A^+1) finite block. Then, 
from the standard recipe [41j to compute the inverse of a matrix!^ another, novel 
expression of the canonical partition function follows: 

ZN^det{nN){Zi)^/N\ (17) 

where f2 jv is an almost triangularized NxN matrix of which the entries are defined 

by 

Aa-b+i/Ai for h<a 

^N\aM-^\ -a/\i for b = a + \ (18) 

otherwise . 



^From det(/ — A) = l, an intermediate relation between (|16p and p7p follows: 

(J- A)-^[iV+l,l] = (-l)^det($iv), 

where $]v is an NxN matrix whose entry $iv[a, b] is given by —Xa-b+i/a for b < a, unity 
for b = a+1 and zero otherwise. 
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In particular, every diagonal entry is unity so that when Ai oo, we have det(r2Ar) — > 
1 and hence the reduction: (Zi)^ /Nl as in ([T^ . 

In addition to the physical quantities 1^ above, by considering —f5~^dE^ InZ = 
■qe"^^^ / [1 — r]e~^-^") as a trick [42], we can also compute the number of particles 
occupying the ground state: 

(^o) = Ek=i e-^P^-Z^^u/ZN ■ (19) 

Each term in the sum above corresponds to the probability for at least k par- 
ticles to occupy the lowest state. If we denote this probability by pk, the dif- 
ference pk — Pk+i corresponds to the probability for precisely k particles to oc- 
cupy the ground state [121111] ■ This leads to an alternative derivation of as: 
(^o) = LfcLi f^iPk - Pk+i) = Y.k=iPk where pn+i = 0. 

Henceforth, exclusively for ideal Bose gas we focus on N particles with mass 
TO, confined in a box of dimension d and length L = V^^'^. Hard, impenetrable 
walls impose Dirichlet boundary condition |45j . Since we are interested in a finite 
system, the periodic boundary condition which is somewhat more popular in the 
literature is not suitable for our purpose. We recall that nevertheless enforcing 
periodic or Neumann boundary condition leads to a thermodynamic instability at 
low temperature near absolute zero for arbitrary N (see Figure 3 in Appendix). 

With positive integer valued good quantum numbers: 

n = (ni,ri2, • • • ,rid) , (20) 
the single particle Boltzmann factor in ^ assumes the form: 

In terms of a Jacobi theta function: 

m---M0,q)-l]/2 = j:Z,q'^\ (22) 
we get specifically for pO)): 

Xk{q) = mq')f. (23) 
After all, Zn becomes a function of only one variable q, so that we may put 

I3dfi = q\nqdq, Vdy = -{2/ d)q\nqdq . (24) 
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This implies that ah the dimensionless physical quantities such as Cy /k^, Cp/k^, 
{Nq), etc. are also functions of the single variable q. Consequently, it turns out 
that the temperature dependence of all these dimensionless quantities can be best 
analyzed if we introduce the following two dimensionless "temperatures" : 

TV k^T{V/N)y''{2m/ {tt^H^)) = (4/7r)(U./A,.B„,He)2 , ^^^^ 
TP k^TP-y^^+^') (2m/(7r2n2))'i/('i+2) ^ 



where Aj^Brogiio — h^J 2t: / {mkB,T) is the thermal de Broglie wavelength and l^ff — 
{ylNY/"^ is the average interparticle distance. While the former in ([25]) is a mono- 
tonically increasing function of the latter may be not so as: 

-\/tv ^ N^/'^nq, -\/Tp = \xiq[{2/d)qdq\xiZ]sf'^'^+'^K (26) 

Any critical value of these quantities will automatically give us the critical temper- 
ature at an arbitrarily given volume or pressure. 

In a similar fashion, we also define a dimensionless indicator of the thermody- 
namic instability, 0, as well as dimensionless "volume" and "energy" : 

^ := -{l/N)liV^dvP = 4C^/[d^k^{Cp - Cv)] , 

VP := (tv/tp)^/^ = (V/N) {2mP/{n^h^)Y'^''^^^ , (27) 

ep := (S/Ar)F-2/('i+2) {2m/ {^'^ti')Y'''''^'''^ . 
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• Low temperature limit: 

The variable q lies between zero and one. As q — > we have 

Zj, ^ e-^^^° = . (28) 
Hence at q = 0, the temperatures vanish ry = rp = 0, and 

Cv = Cp^O, (/) = oo , {No)^N, vp^ (2d/Ar2)i/(<i+2) _ (29) 
In particular, the volume reads at absolute zero, 

V = [NTT^h^l{mP)Y''^'^+'^^ . (30) 

That is to say, despite of the apparent Bose-Einstein condensation i.e. (Nq) = N, 
the volume assumes a finite value which is even not extensive. The finiteness is 
essentially due to the Heisenberg uncertainty principle: Since the particles are lo- 
calized in a finite box, the uncertainty principle forbids the ground state energy Eq 
to vanish and leads to the nontrivial canonical partition function (|28p . 

• High temperature limit: 

As q —?■ 1, from (|17p and thanks to an identity of the theta function by Jacobi 

^?(e-"") + 1/2 = [i?(e-"/") + l/2]/V^, (31) 

we have 

Zn (l/Af!)(-7r/41ng)^^/' . (32) 
Hence at q = 1, the temperatures diverge Ty = Tp = oo, and 

Cv/kn^d/2, Cp/k^ = l + d/2, 0=1, 

(33) 

vp/rp = 1 , {Tvr/(^+^^/Tp = 1 , (iVo) = . 

Namely, ideal Bose gas reduces to the classical ideal gas at high temperature, sat- 
isfying the classic relation: 

PV = Nk^T . (34) 
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3.2 Numerical results 



Here we present our numerical results of o? = 3, i.e. cubic box, at generic temper- 
ature. Our analysis is based on a set of recurrence relations which enables us to 
perform iV^ order computation, enhanced by a parallel computing power - JS20 
(PowerPC 970) system. Essentially we utilize p"3)FI and two other relations which 
can be straightforwardly obtained from the grand canonical partition function ex- 
pressed in the form (flUl) : 

N 

PdpZj^ — ^ Pn\nZN-n ) 

(35) 

P'^d'pZN = ^ [n{Cn + p'i)\nZN-7i + PnK Pd^ZN-n] , 

where we set 

/9„ (/3/n)9^ In A„ , (« := (/3V"')9^ In A„ . (36) 

All our results agree with the asymptotic behaviours ((29|) . ((33)) . 

• Constant volume curves (Figure 1): 

As expected, all the physical quantities are smooth single valued functions of the 
temperature ry. As grows, the specific heat Cy develops a maximum, at which 
(A^o) drops rapidly. This behaviour is consistent with Ref. [24]. More importantly 
for us, while (/) decreases from infinity at Ty = to one at ry = oo, it develops a 
local minimum which becomes eventually negative if N > 7616. This manifests the 
thermodynamic instability of the ideal Bose gas confined in a cubic box. 



^Since Zn can be a big number for large A'^, for the evaluation of the canonical partition 
function itself, it is convenient to decompose it as 

= nr=i ' 

and utilize a recurrence relation: 

/iV=E:roMA^-n/(iVnf_V,/,)], 

which is equivalent to (|13p . 
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Cv/kB <No>/N 




Figure 1: Constant volume curves 

a: specific heat per particle at constant volume, Cy/k^ 

b: the occupancy ratio of the ground state, {Nq)/N 

c, d: the thermodynamic instability indicator, cj) = —{l/N)l3V^dYP 

Here a,b,c are for A^= 1 (red), iU (orange), 10"^ (green), 10"^ (blue), while d 

is for = 7614, 7615, 7616, 7617. The case of N= 1 (red) also corresponds 

to ideal Boltzmann gas. All the quantities are dimensionless. 
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Figure 2: Constant pressure curves 

a: the dimensionless volume vp versus the dimensionless temperature rp 

b: Ty versus Tp 

c: {Nq)/N versus Tp 

d: the dimensionless energy ep versus rp 

e, ee: specific heat per particle under constant pressure Cp/k^ versus Tp 
Note that a,b,c,d,e are for A^= 1 (red) i.e. ideal Boltzmann gas or 
N= 10 (orange), 10~ (green), 10^ (blue). The small boxes magnify the zigzag 
segments of = 10^ (blue), while ee magnifies e for = 10^ (blue). The 
dotted pink lines denote supercooling and superheating points. 
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• Constant pressure curves (Figure 2): 

Since the thermodynamic instabihty indicator (j) vanishes at two points when N > 
7616, there are genericahy two critical temperatures, Tp < Tp* , which we suggest 
to identiiy as supercoohng and superheating points respectively. As predicted from 
the general argument ([5]), all the physical quantities change discontinuously at 
these points: On each physical quantity versus temperature plane, the constant 
pressure line zigzags between the two critical points keeping them as two turning 
points. Accordingly, physical quantities are triple valued between the two points 
and double valued at the points. This implies the existence of three different phases 
during the phase transition, say two 'forward' phases and one 'backward' phase. 

Specifically, the volume expands abruptly (Figure 2a): below the supercooling 
point it is almost constant whilst beyond the superheating point it follows the 
classical ideal gas law which resembles the usual liquid-gas transition. At the 
same time, the number of particles on the ground state drops abruptly from the 
full occupancy to the total evacuation (Figure 2c). This clearly realizes a Bose- 
Einstein condensation taking place both in the momentum and in the coordinate 
space [47j. The specific heat per particle under constant pressure diverges at the 
critical points, and when magnified between the supercooling and the superheating 
points (Figure 2ee), intriguingly it reveals one negative heat capacity in addition 
to other two positive ones. Because Cp/k^ — dep/drp + d-up/drp, discontinuous 
changes in both the internal energy and the volume simultaneously contribute to 
the divergence of the specific heato 

Algebraically, while physical quantities are a priori single valued functions of 
q and hence ry, the origin of the multi-valuedness and singular behaviours can 
be all traced back to the zigzagging of the isobar curve on (Ty,rp) plane (Figure 
2b). Namely the isobar curves of the ideal Bose gas zigzag on the (V,T) plane, if 
N > 7616. 

Our numerical results of the supercooling and the superheating points are listed 
below for selected N: 

*In fact, the exponent corresponding to the singularities can be shown to be 1/2 [48j . 
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N 


T* 


** 


7616 


1.0543694113 


1.0543694116 


10* 


1.05270 


1.05277 


10^ 


1.0410 


1.0424 


lO^ 


1.034 


1.036 



from which, recovering all the dimensionful parameters, we may read off the super- 
cooling and the superheating temperatures at arbitrary constant pressure for each 
value of N: 

h T ~ T* V r ^"^^ ^13/5 02/5 

'i'B supercool — ' P \ Im ) ' 

(37) 

h T T** V ('2£ft^^3/5p2/5 

"-B superheat — 'P V ~2m~ ' 

With the mass of Helium-4 and under the pressure 1 atm, our numerical result 
of N = 10^ gives us the supercooling temperature 1.686 Kelvin and the superheat- 
ing temperature 1.689 Kelvin, both of which are of the same order as the lambda 
point 2.17 Kelvin or alternatively as the boiling point 4.22 Kelvin. 



4 Discussion 

In summary, finite amount of ideal Bose gas confined in a cubic box reveals a ther- 
modynamic instability, if TV > 7616. This implies that, binder constant pressure 
condition, the ideal Bose gas may undergo a first-order phase transition accom- 
panied by genuine mathematical singularities. It is characterized by both the su- 
percooling and the superheating points; Bose-Einstein condensation in both the 
momentum and the coordinate spaces; and the very fact that all the physical quan- 
tities become triple valued between the two points. In particular, one of the three 
values of the specific heat under constant pressure is negative. While Bosc-Einstcin 
condensation appears as a continuous phase transition when the volume is kept 
constant (Figure lb), it may become a discrete phase transition if the pressure is 
held fixed (Figure 2 c). 

Ideal gas consists of featureless point particles. It is such a simple model that 
the canonical partition function becomes essentially a one- variable function, Zi^{q). 
Assuming an extra internal structure of the particle, e.g. spin or a vibrational mode. 
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will introduce an additional dimensionful parameter into the energy spectrum. This 
may generate an additional (continuous or discontinuous) phase transition and dis- 
tinguish the Bose-Einstein condensation to the ground state from the condensation 
in the coordinate space [35] . 

The critical number 7616 we report in this paper can be regarded as a char- 
acteristic number of 'cube' that is the geometric shape of the box containing the 
ideal Bose gas. Boxes of different shapes (see e.g. BQ ) will have different critical 
numbers. Thus, our scheme of investigating the thermodynamic instability of ideal 
Bose gas can provide a novel algorithm to assign a characteristic number to each 
geometric closed two-dimensional manifold. 

Apparently ideal Bose gas has no interaction. However, compared to ideal 
Boltzmann gas, ideal Bose gas has higher probability for the particles to occupy 
the same quantum state, as seen from the comparison between (HH) and (HUli In 
coordinate space this means that identical bosonic particles have tendency to gather 
together in comparison to ideal Boltzmann gas. In fact, a path integral representa- 
tion of the canonical partition function of ideal Bose gas reveals an existence of an 
attractive effective potential [6 |f| This provides a physical clue to the condensation 
in both the momentum and the coordinate spaces: As the temperature decreases, 
the effective statistical attraction becomes dominant and the system condensates. 
In this context, it is also worth while to recall the similarity between the permu- 
tation symmetry of the identical particle indices and the gauge symmetry of the 
standard model in high energy physics or matrix models (see e.g. \52\). Although 
the former is discrete while the latter is continuous, the latter may include the 
former as a subgroup. The common feature is that they both correspond to non- 
physical symmetry [53 . As a matter of fact, up to an overall factor, the canonical 
partition function of identical harmonic oscillators [37] coincides with that of a 
massive Yang-Mills quantum mechanics [54], taking the form n^^=i(l ~ 

Generically, for a stable matter dyP is negative. What we show by taking ideal 

^As a simple example, consider a two-particle system with quantum states 'up' and 
'down' only. The probability for the two identical bosonic particles to occupy the same 
'up' state is 1/3 and this is higher than that for the case of two Boltzmann particles, 
1/4. Similarly, the probability for the two identical bosonic particles to occupy the two 

different states is 1/3 which is lower than that of the two Boltzmann particles, 1/2. 

®It is worth while to note that, the nontrivial volume at absolute zero due to the 

Heisenberg uncertainty principle (|30|) suggests a statistical repulsive interaction at low 
temperature. Nevertheless this is valid for ideal Boltzmann gas too. For related subtle 
issues see e.g. [51| . 
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Bose gas as an exactly solvable model is an explicit demonstration that, if there 
are sufficiently, yet finitely, many identical bosonic particles, dyP can be positive 
for at least one interval of temperature. It will be therefore interesting and crucial 
to see, to what extent intciractions can alter this. If not much, as in a weakly 
interacting system, one first-order phase transition accompanying a discontinuous 
volume change, or the liquid-gas transition itself, is likely to occur essentially due 
to the identical nature of particles. In a perturbative analysis of the interaction, 
ideal Bose gas provides the 'zeroth' order contribution: For each quantum state tp, 
as the energy eigenvalue is shifted from to + A^, the canonical partition 
function changes from Zjv to Zj^{e~^^'i') . Thus, the thermodynamic instability is 
modified as 

cj) ^ cj)- {l/N)V^d^ylu{e-^'^^) = <^ + {l/N)pV^dl{I^^) +■■■ . 

If the first order correction gives a negative contribution, 9y(A^) < 0, the ther- 
modynamic instability persists naturally up to the order. In this case, the average 
bonding energy between particles will also change discontinuously during the phase 
transition, but this will be an accompanying side effect rather than a key reason to 
drive the phase transition. 

It will be experimentally challenging to find a corresponding critical number for 
each molecule to manifest a discontinuous phase transition or its liquid-gas transi- 
tion under constant pressure. A criterion for the first-order phase transition is to 
observe the supercooling and the superheating phenomena. 
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Appendix: Other boundary conditions 



Here as an explicit demonstration of the theorem (|8]), we present numerical results 
(Figure 3) on the ideal Boltzmann gas (or equivalently ideal Bose gas with N= 1) 
in a cubic box under boundary condition: Neumann or periodic. This corresponds 
to taking Afc(g) = [1 + 'd{q'')]'^ for the former and Xk{q) = [1 + 2i){q'^'')]'^ for the 
latter. Actually we only need Xi{q) and set d = 3 as for the cubic box. 

As expected, due to the existing zero mode, both boundary conditions lead to 
a thermodynamic instability i.e. dyP > at low temperature near absolute zero. 




Figure 3: Thermodynamic instability due to zero mode 

a: (j) versus Ty under Neumann boundary condition 
b: (p versus ry under periodic boundary condition 
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